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The mixing of photons with axion-hke particles ( ALPs) in the large-scale magnetic field B changes 
the polarization angle of a hnearly polarized photon beam from active galactic nuclei in radio galaxies 
as it propagates over cosmological distances. Using available ultraviolet polarization data concerning 
these sources we derive a new bound on the product of the photon- ALP coupling ga-y times B. We 
find ga-yB < 10"" GeV"^ nG for ultralight ALPs with m^ < 10"^^ eV. We compare our new 
bound with the ones present in the literature and we comment about possible improvements with 
observations of more sources. 
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I. INTRODUCTION 

Many extensions of the Standard Model predict the ex- 
istence of axion-like particles (ALPs) , namely very light 
spin-zero bosons a characterized by a two-photon ver- 
tex a77 T]. In the presence of an external electromag- 
netic field, this coupling gives rise to photon- ALP mixing, 
which leads in turn to two distinct phenomena. One is 
photon- ALP oscillations [2H1], which is exploited by the 
ADMX experiment [5] to search for ALP dark matter [5] , 
by CAST to search for solar ALPs [7] and by the regen- 
eration laser experiments [8lfT2] . The other phenomenon 
consists in the change of the polarization state of a photon 
beam traveling in a magnetic field. More specifically, an 
initially linearly polarized beam propagating in a mag- 
netic field undergoes two different effects: (a) it acquires 
an elliptical polarization, and (b) the ellipse's major axis 
gets rotated with respect to the initial polarization direc- 
tion. The former effect is called birefringence while the 
latter dichroism [3l [131 HI] • A claim for a positive obser- 
vational evidence of both vacuum birefringence and vac- 
uum dichroism by the PVLAS collaboration [TS] employ- 
ing a laser beam has subsequently been withdrawn [16] . 

Over the last few years it has been realized that the 
077 coupling can also produce detectable effects in the 
observations of distant active galactic nuclei (AGN) and 
gamma-ray bursts (GRBs), since photons emitted by 
these sources can mix with ALPs during their propa- 
gation through large-scale magnetic fields (see T] for 
a recent review). In this context, photon- ALP oscilla- 
tions |17H20) provide a natural mechanism to drastically 
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reduce photon absorption caused by extragalactic back- 
ground light above roughly 100 GeV where observations 
of AGN are currently performed by Imaging Atmospheric 
Cherenkov Telescopes 21 29j. Astrophysical observa- 
tions of AGN and GRBs also offer the opportunity to 
look for ALP effects in polarization measurements. Con- 
straints on ALPs from future GRB polarimctric measure- 
ments have been recently discussed [301 IS] • 

In this paper we derive a new bound on the photon- 
ALP coupling using ultraviolet (UV) data from AGN in 
radio galaxies (RGs) by searching for a rotation of the 
detected photon linear polarization with respect to the 
emitted one. Our analysis has been inspired by recent 
studies [311 [33] , where the same data-set has been used 
to constrain the effects of a possible cosmological dichro- 
ism, occurring when photons coupled to an hypothetical 
background of cosmic Nambu-Goldstone fields propagate 
over long distances. However, our investigation is quite 
different in spirit. For, we do not absolutely assume the 
presence of a background of ALPs but merely the exis- 
tence of a large-scale magnetic field B. 

Now, in spite of the fact that AGN in RGs are very 
complicated objects, the scattering model of anisotropic 
nuclear UV radiation predicts that the angle between the 
UV axis of the source and the UV linear polarization is 
90° [32]. Cosmic dichroism or photon- ALP mixing can 
- by rotating the linear UV polarization - give rise to 
an angle different from 90° between the observed linear 
polarization of the UV radiation and the UV axis of a sin- 
gle source. Given that no rotation has been detected so 
far within a few degrees for each analyzed source, we are 
able to set a bound on the product of the photon- ALP 
coupling times the strength of the large-scale magnetic 
field - namely on ga-yB - since any observable effect in- 
volving ALPs in an external magnetic field depends on 
the product ga-yB. 

The plan of the paper is as follows. In Sec. 2 we re- 



view the formalism describing the photon-ALP mixing 
and its implications for photon oscillations and polariza- 
tion. In Sec. 3 we introduce our input parameters, the 
AGN in RGs data-set and we describe the constraints 
coining from the absence of rotation of the linear polar- 
ization of the UV radiation emitted by the considered 
sources. Finally, in Sec. 4 we compare our new bound 
with the the ones present in the literature, and we dis- 
cuss possible improvements. 



II. MIXING OF PHOTONS WITH AXION-LIKE 
PARTICLES 

A. Oscillations 

Photon-ALP mixing occurs in the presence of an ex- 
ternal magnetic field B due to the interaction term [2]-[4] 
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where ga-y is the photon-ALP coupling constant (which 
has the dimension of an inverse energy) . 

We consider throughout a monochromatic and linearly 
polarized photon beam of energy E propagating along the 
z direction in a cold ionized and magnetized medium. It 
has been shown that for very relativistic ALPs, the beam 
propagation equations takes the form [3j 
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where Ax{z) and Ay{z) are the photon linear polarization 
amplitudes along the x and y axis, respectively, a{z) de- 
notes the ALP amplitude and Ai represents the photon- 
ALP mixing matrix. 

Actually, A4 takes a simpler form if we restrict our 
attention to the case in which B is homogeneous. We 
denote by By the transverse magnetic field, namely its 
component in the plane normal to the beam direction. 
We can choose the y axis along By so that B^ vanishes. 
The linear photon polarization state parallel to the trans- 
verse field direction By is then denoted by A\\ and the 
orthogonal one by A±. In this case the mixing matrix 
can be written as fSH 551 
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whose elements are i3J Aj_ = Api -I- A^^, A|| = Api + 
ACM, ^^^ = 1/2.9,^St, A„ = -ml/2E. The term 
A.- = -'-2 
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is the plasma frequency expressed as a function of the 
electron density in the medium Ue. The terms AP^ de- 
scribe the Cotton-Mouton effect, i.e. the birefringence 
of fluids in the presence of a transverse magnetic field 
where lA?^ - A9^\ cx BI. The latter terms are of little 



importance for the following analysis and will thus be ne- 
glected. The Faraday rotation term Ar, which depends 
on the energy and the longitudinal component Bz, cou- 
ples the modes A\\ and A± [36] . While Faraday rotation 
is in principle important when analyzing polarized pho- 
ton sources, it actually plays no role at the energies we 
are dealing with. 

In this simplified geometry, the component A± decou- 
ples away, and the propagation equation reduce to a 2- 
dimensional problem. Its solution follows from the di- 
agonalization of the 2-dimensional mixing sub-matrix of 
A^o through a similarity transformation performed with 
an orthogonal matrix, parametrized by the rotation angle 
6 which takes the value |3] 
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In particular, the probability for a photon emitted in 
the state Ay to oscillate into an ALP after traveling a 
distance d is given by [3] 



r> • 2o/l • 2 / ^osc d 

Pa-y = sm 20 sm ' 



(Aajd) 



sin2(Ao,ed/2) 



where the oscillation wave number is 
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Aosc= [(A,-Api)' + 4A2^] 
It proves useful to define a critical energy 
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in terms of which the oscillation wave number can be 
rewritten as 




(9) 



From Eqs. (l5])-((9]) it follows that in the energy range E ^ 
Ec the photon-ALP mixing is maximal {9 ~ 7''/4) and 
the conversion probability becomes energy-independent. 
This is the so-called strong-mixing regime. Outside this 
regime the conversion probability turns out to be energy- 
dependent and vanishingly small, so that Ec acquires the 
meaning of a low-energy cut-off. As already pointed out, 
it is now evident that all results depend on the product 
gajB, and not on ga-y and B separately. 



B. Polarization 

Let us now focus our attention on the polarization ef- 
fects, keeping in mind that for our monochromatic beam 
the electric field components are £^11 oc Au and E± ex A±. 
We also denote by 7|| and 7_l the beam photons described 
by A|| and A±, respectively, and by jb a "magnetic field 
photon". Since the coupling in Eq. (IT]) has the form 
A|| B-T a, it follows that the photons 7|| mix with the 
ALPs a whereas the photons 7_l do not mix at all. 

Assuming the considered beam to be initially {z ~ 0) 
linearly polarized at an angle ipo (with respect to an arbi- 
trarily chosen fixed direction), we address first the effect 
arising from the production of real ALPs, occurring via 
the 077 vertex with one photon line representing a 75 
photon. Since only 7|| mixes with a, the || mode gets de- 
pleted thereby decreasing Ay whereas A± is manifestly 
unchanged. As a consequence, we have a rotation of A as 
the beam propagates, and since it occurs because of the 
selective loss of photons depending on their polarization 
the effect in question is called dichroism, owing to the 
analogous situation occurring in classical optics. Clearly 
after a distance z the rotation angle is 



ip{z) = ipo + arctan 



A^{z) 



(10) 



Consider next what happens for the exchange of vir- 
tual ALPs, namely for the diagram in which two 077 
vertices are joined together by the a line and two exter- 
nal photon lines are actually 75 photons. Evidently this 
diagram is nonvanishing for incoming 7|| photons while 
it vanishes for incoming 7_l ones. Consequently, the in- 
dex of refraction in the two modes || and _L is different, 
which means that the two modes 7|| and 7_l travel at 
different speeds. Therefore, at any finite distance from 
the source, the beam polarization turns out to be ellipti- 
cal. Even more explicitly, as the photon beam propagates 
its A-vector changes both direction and magnitude so as 
to trace a spiral around the z direction with elliptical 
sections. After each 27r rotation, a different projected el- 
lipse gets singled out in the plane perpendicular to the 
beam. This effect is called birefringence due to the similar 
situation taking place when a linearly polarized photon 
beam traverses an anisotropic medium. Still - as long as 
birefringence alone is concerned - all such ellipses have 
parallel major axes, which is just an elementary manifes- 
tation of the composition of two harmonic motions along 
orthogonal directions. 

So, when both dichroism and birefringence are at work 
an initially linearly polarized photon beam acquires an 
elliptical polarization with the ellipse's major axis ro- 
tated with respect to the initial polarization. We remark 
that a qualitatively similar situation occurs in vacuum 
according to quantum electrodynamics, where birefrin- 
gence arises from Delbriick scattering (photon scattering 
in a magnetic field) |37j while dichroism is produced by 
photon splitting ^3H] (see |3S] for a review). 



Elsewhere the birefringence in the emission of GRBs 
has been analyzed }30j . whereas here our interest is fo- 
cussed on dichroism. In spite of the difference in the 
physical effects, both studies are based on the propaga- 
tion of photons from distant sources through large-scale 
magnetic fields. According to the standard lore, they 
are modeled as a network of magnetic domains, with size 
equal to the coherence length. In every domain the mag- 
netic field B is assumed to have the same strength but 
its direction is allowed to change randomly from one do- 
main to another. Therefore the propagation over many 
magnetic domains is a truly S-dimcnsional problem, be- 
cause - due to the randomness of the direction of B - the 
mixing matrix A4 entering the beam propagation equa- 
tion cannot be reduced to a block-diagonal form similar 
to A^o in all domains. Rather, we take the x, y, z coor- 
dinate system as fixed once and for all, and - denoting 
by tpn the angle between Bt and the y axis in the n-th 
domain - we treat every ^„ as a random variable in the 
range < ip < 2TT. The numerical technique adopted to 
solve the beam propagation equation in the present case 
closely follows the one used in |3nj , to which the reader is 
addressed for more details. We stress that because of the 
random orientation of the magnetic field in each domain, 
the effect of photon-ALP mixing on the photon polar- 
ization strongly depends on the orientation of the line of 
sight. As a consequence of the stochastic nature of this 
process, the photon polarization will be characterized by 
a probability distribution function, obtained by consid- 
ering photon-ALP conversions over different realizations 
of the large-scale magnetic fields. 



III. AGN CONSTRAINTS 

A. Input parameters 

As a preliminary step, it is convenient to express the 
relevant physical quantities in units set by their nat- 
ural values. The energy of UV photons is ^ 10 eV. 
The strength of the large-scale magnetic field in the 
intergalactic medium has to meet the constraint B < 
2.8 X 10~'^(L/Mpc)"i/2G - where L denotes its coher- 
ence length - which arises by scaling the original bound 
from the Faraday effect of distant radio sources [301 SI] 
to the now much better known baryon density mea- 
sured by the Wilkinson Microwave Anisotropy Probe 
(WMAP) mission [42]. Its coherence length is expected 
to lie in the range 1 Mpc < L < lOMpc JT- . The 
mean diffuse intergalactic plasma density is bounded by 
rig ^ 2.7 X 10^^ cm~'^, arising from the WMAP measure- 
ment of the baryon density [42]. This corresponds to a 
plasma frequency LUp\ < 1.8 x 10""'^'* eV today. However, 
at the moment no measurements of n^ are available in 
the large voids of the interstellar medium. Since the av- 
erage value of Ue inferred by WMAP has been sometimes 
criticized [H] claiming it to be too large for most of the 
intergalactic space, we allow it to vary within a range of 



two orders of magnitude below this value. Accordingly 
the upper bound on Wpi gets reduced by one order of mag- 
nitude. Recent results from the CAST experiment yield 
a direct bound on the photon-ALP coupling constant 
ga-y < 8.8 X 10-" GeV-i for m < 0.02 eV iZj, slightly 
better than the long-standing globular-cluster limit [45] . 
In addition, for m < 10"^° eV a more stringent limit 
arises from the absence of 7-rays from SN 1987A, giv- 
ing ffa7 < 1 X 10"" GeV"^ gnillTj even if with a large 
uncertainty. Therefore, a suitable parametrization of the 
relevant quantities is 
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Moreover, the value of the critical energy in Eq. (|8| is 
given by 
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As we are interested in studying the propagation of 
photons from AGN at redshifts z ~ 3 we must also take 
into account the redshift evolution of the environment 
that causes the photon-ALP mixing. Notably the ob- 
served energy Eq is red-shifted with respect to emit- 
ted one E, i.e. E — Eq{1 + z). Assuming that the 
large-scale magnetic field is frozen in the plasma [H], 
its strength scales as B{z) = Bo(l + z^ and clearly we 
have ne{z) — nefi{l + zY . Here, the subscript indicates 
values at the present epoch. The physical size of a mag- 
netic domain scales as L{z) — Lq{1 -^ z)^^ provided that 
it is smaller than the Hubble radius. 



TABLE I. Linear UV polarization in distant AGN in RGs [32]. 



Name 


Redshift {z) 


ARA. (deg) 


MRC 0211-122 


2.34 


89.0 + 3.5 


4C -00.54 


2.363 


82 + 8 


4C 23.56a 


2.482 


94.6 + 9.7 


TXS 0828+193 


2.572 


91.6 + 4.5 


MRC 2025-218 


2.63 


86 + 9 


TXS 0943-242 


2.923 


89.7 ±4.4 


TXS 0119+130 


3.516 


95 + 16 


TXS 1243+036 


3.570 


86.0 + 8.8 



the orientation of the polarization of optical/UV radia- 
tion can be predicted from physical arguments [3H [15] 
and it is expected to be perpendicular to the optical/UV 
source axis. Actually, astronomers define a fiducial fixed 
direction on the plane of the sky and define the source 
position angle as that formed by the optical/UV source 
axis with such a fiducial direction. Since photon polar- 
ization is transverse, the observed ones have polarization 
direction in the plane of the sky, and so one can define 
the polarization position angle as the one formed by the 
polarization and the above fiducial direction. So, we can 
rephrase the above expectation by stating that the dif- 
ference AP.A. between the source and the polarization 
position angles is expected to be 90°. 

We consider here and show in Table[l]the set of sources 
already used in [32] to constrain the cosmological dichro- 
ism. The second column indicates the redshift z of each 
source. The third column reports AP.A., which has been 
measured on the available images. These sources were 
selected in order to provide robust upper limits on the 
rotation of the polarization angle in the UV band (1300 
A in the source rest-frame) [35] ■ Indeed, the fact that 
AP.A. is close to 90° for every object - actually compat- 
ible with 90° within the accuracy of the measurements 
- puts stringent constraints on the rotation of the linear 
polarization of the UV radiation during the propagation. 



B. Data 



In order to quantify the amount of rotation of the linear 
polarization which takes place during the propagation, it 
is important to deal with photon sources for which the 
initial photon linear polarization can be reliably deter- 
mined. Optical/UV radiation from AGN in RGs at red- 
shift z > 2 is very well suited in this respect. These 
astronomical objects are well known emitters of linearly 
polarized radiation in the radio band as well as in the 
optical/UV bands. While the dichroism analysis of radio 
photons suffers from strong uncertainties, coming both 
from the absence of a reliable prediction for the polariza- 
tion characteristics and from Faraday rotation, the analy- 
sis of optical/UV radiation is more robust. In particular. 



C. Likelihood analysis 

We derive our constraint on the rotation of the polar- 
ization angle by applying a likelihood technique accord- 
ing to a standard procedure }49j . As it is clear from the 
discussion presented above, this effect can depend only 
on two independent quantities: gajB and n^, which then 
span our parameter space. In order to reconstruct the 
probability distribution function of the "signal" in the 
presence of photon-ALP mixing we perform 10^ Monte 
Carlo runs of the beam propagation from the source to 
us in the intergalactic medium - with the properties con- 
sidered above - for a given point in the parameter space 
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FIG. 1. Probability density functions p{Aip) plotted versus 
the rotation of the polarization angle Aip after propagation 
in the extragalactic magnetic field, simulating 10^ AGNs at 
redshift z = 3 with initial polarization angle ipo = 90° . We ad- 
dress the case with z-evolution and we take n^ = 10~* cm~'^. 
Three representative values of Qa-yB are considered. 



{ga-yB, Ue). We take z = 3 for the redshift of the source, 
which corresponds to the typical redshift of the sources in 
Table |lj In each simulation, the large-scale magnetic field 
is modeled as a network of random magnetic domains. 
We consider both the situation in which the magnetic 
field and the electron density evolve with the redshift (z- 
evolution case) and the one in which no redshift evolution 
takes place (no z-evolution case). We compute then the 
likelihood of the hypothesis of photon-ALP mixing by 
evaluating 



^ALP ^^^=lPi^^,(^^\Q) , 



(15) 



where ipi, i — 1...N denote the data on the rotation of the 
polarization, ai are the errors associated with the indi- 
vidual measurements and = (ga-yB^rie) represents the 
parameter space to be constrained. The function p{x\y) 
is the probability distribution function obtained with the 
Monte Carlo simulations, weighted with the experimen- 
tal resolution (assumed Gaussian) associated with the 
measurement ipi 
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In Fig. [T] we show the dependence of the probability 
density functions p{Aip) on the rotation of the polariza- 
tion angle Aip = ip{z) — tpo with respect to the initial 
value ipo = 90°, assuming for simplicity perfect resolu- 
tion R — 1. As an illustration, we address the case with 
z-evolution and we take n^ — 10^* cm^^. In the Figure, 
for ga-fB — 10"^"^ GeV~^ nG the effect on the polar- 
ization is negligible, the probability distribution being 
a delta-function peaked on Aip = 0. For a larger val- 
ues of ga-yB = 3 X 10~^^ GeV~^ nG, the mixing tends 
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FIG. 2. Exclusion plot for photon-ALP conversion based on 
UV polarization data from AGN in RGs. The colored regions 
are excluded at 95 % CL. In particular, the region above the 
continuous curve is excluded for 2-dependent case, while the 
region above the dashed curve panel is excluded for the z- 
independent case. 



to broaden the AGN polarization distribution around 
A^ = 0. Finally, for ga-yB = 10"!° GeV"! nG, the 
probability distribution is completely flat in Aip. 

The likelihood for the null hypothesis (no mixing) 
reads instead 



inuU = nj^;^i?((y9j,(Ti) 



(17) 



Assuming that Wilks' theorem applies [50], the confi- 
dence region compatible with the null rotation hypothesis 
is given by the condition 



IuLalp - hiLnuu < -Xp{k) 



(18) 



where k is the number of parameters, i.e. fc = 2. Note 
that Xfs{k) = 2.3, 4.61 and 6.17 for /? = 68.3%, 90%, and 
95.4%, respectively. 

In Fig. [2] we show our exclusion contour in the plane 
(ng, gajB). The colored regions are excluded at 95 % 
CL. In particular, the region above the continuous curve 
is excluded in the z-evolution case, while the region above 
the dashed curve is excluded in the no z-evolution case. 
In the former case the values of Ue and ga-yB shown on 
the axes are to be understood at z = 0. Within a factor 
of a few, the same contours also hold if one increases the 
domain size by a factor of 10. 

We find ga^B < 2 x 10"" GeV"! nG in the no z- 
evolution case, and the bound improves in the z-evolution 
case reaching ga-yB < 5 x 10~^^ GeV~^ nG for n^ < 
7 X 10"^ cm"3. 
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FIG. 3. Bounds for the product of ALP coupling to photon 
times magnetic field Qa-yB as a function of the average electron 
density rie. The gray regions are excluded from the SNIa 
dimming for the jz-independent case (dark) and 2-dependent 
case (light) [52]. Also shown are the regions constrained by 
CMB [34] (above the black dashed line), which dominate at 
low Ue and QSO [53) (light blue) spectral distortions. The 
UV exclusion region is the one above the solid curve for the 
z-dependent case, and above the dashed horizontal line for 
the z-independent case. 



IV. DISCUSSION AND CONCLUSIONS 

We have found a new limit on the quantity ga-jB for a 
ultralight ALPs with rua < 10^^^ eV by studying the ro- 
tation of the linear polarization measured in UV photons 
from AGN in radio galaxies at cosmological distances. 

It is now interesting to compare our bound with other 
ones recently discussed in the literature, concerning ul- 
tralight ALPs and based on photon-ALP oscillations in 
the photon beam from distant sources in large-scale mag- 
netic fields. We show in Fig. [3] a summary of these limits. 

The photon-ALP oscillations can lead to a dimming 
of the Type la supernovae luminosity curves [HI IFT] . 
This effect has been recently studied to place bounds 
on QajB [55]. Our limit is stronger by roughly a fac- 
tor of five than the one resulting from this effect for 
Tie ^ 10""^ cm~'^, while for Ue ^ 10^* cm""^ the two 
bounds become comparable in the no z-evolution case 
but our one is slightly better in the z-evolution case. 

Photon-ALP oscillations also imprint peculiar distor- 
tions on quasar (QSO) colors and spectra [53]. The cor- 
responding bound obtained in the no z-evolution case is 
worse than our by a factor of three for rig > 10"^ cm^'^ 
and becomes slightly better for lower values of the elec- 
tron density. 

Finally, photon-ALP conversions can cause an exces- 
sive spectral distortion of the cosmic microwave back- 
ground (CMB) |34l |54j- However, this limit would be 



relevant only in the case of extremely low electron den- 
sity, i.e. smaller than ^ 10~^" cm~'^, giving ga-yB < 
fewx 10^^^ GeV^^ nG, or when the ALP is not ultralight 
and resonant photon-ALP conversion is possible [51j . 

From our bound on ga^B it would be possible to obtain 
a constraint on the photon-ALP coupling gaj- However, 
due to the very uncertain strength of the large-scale mag- 
netic field B, this bound is uncertain by several orders of 
magnitude. If the large-scale magnetic field had a value 
close to the upper bound _B ~ 10 nG for a domain size 
Lg — 1 ~ lOMpc [55], our astrophysical bound for ul- 
tralight ALPs would improve by at least two orders of 
magnitude the experimental limit from CAST [T^] and 
by at least by one order of magnitude the one obtained 
from SN1987A, g^^ < 10"" GeV"! ^. Recently, ob- 
servational evidence has been discussed that the inter- 
galactic medium is efficiently heated through generation 
of plasma instabilities by powerful blazars [56 . If this 
heating mechanism is at work, it would imply a model- 
dependent upper limit on the the large-scale magnetic 
field of J5 ^ 10"'^ nG. This would imply the bound 
on the photon-ALP coupling g^-y ^ 10"^ GeV~^, sim- 
ilar to the current one from the "light shining through 
the wall" experiment ALPS at DESY, namely ga-y ^ 
7 X 10~® GeV~^ ^2\- Conversely, if the large-scale mag- 
netic field were close to the lower bound B ~ 10~^ nG 
recently inferred from Fermi-LAT data [ 57] . this would 
imply a bound on ga-y < 10^^ GeV^^. It is interesting 
that in this case our limit - based on a small set of data 
and obtained without performing a dedicated study - 
would be comparable with the bound resulting from the 
laser polarization experiment performed by the PVLAS 
collaboration [16] . 

Future improvements of our bound are expected when 
more data on UV polarization from distant AGN in radio 
galaxies will be available. In particular, we have tested 
the effect of improving the precision in the determina- 
tion of rotation and the effect of increasing the number 
of detected sources. For an error of order of 2° in the de- 
termination of absence of rotation (to be compared with 
the ~ 6° error in present data) we find only little im- 
provement. On the other hand, increasing the number 
of detected sources from 8 to 20, and considering again 
an error of ^ 2° would improve by roughly a factor of 
two the constraint on ga-fB. In this respect, we remark 
that the number of the considered sources at z > 2 is 
close to 200 [5S] of which 10 objects have been observed 
with spectro-polarimeters and 8 have been found to be 
polarized. Therefore, it is not unrealistic to expect an im- 
provement in the statistics of the UV polarized sources 
in question in the near future. 

In conclusion, our work suggests once more the inter- 
esting physics potential of astrophysical observations as a 
means to constrain ALPs. From the above discussion it 
follows that all the current bounds on ultralight ALPs 
nicely merge. Nevertheless, since every experimental 
measure and every cosmological and astrophysical con- 
straint has its own systematic uncertainties and its own 



recognized or un-recognized loopholes, it is important to 
use many different approaches in order to constrain these 
elusive particles. 
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